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The time dependent rotation of one orthogonal coordinate system with respect to a fixed one is 
considered in the parametrization based on the effective axis and angle of the rotation, a parametri- 
zation which has recently been used to discuss the irreducible representations of the rotation group. 
The method of the intrinsic vector is used to derive the equations of motion for the instantaneous 
effective axis and angle. A new representation of the angular velocity is obtained in a rotating or- 
thogonal coordinate system generated by a unit vector along the effective axis, and a new geometrical 
interpretation of the effective angle is given. 
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1 . Introduction 

The most common parametrization of the rotation group is by means of the Euler angles 
[1, 2]. 1 Recently, there has been some interest in the parametrization of this group by the direct 
use of the effective axis and angle of a rotation. Moses [3, 4] has calculated the irreducible represen- 
tations in that parametrization, and [5] has computed the orthogonality relations between the 
matrices of any two such irreducible representations. Carmeli [6] has obtained the same results 
using a technique which differs from Moses', and which is originally due to Weyl. The outstand- 
ing feature of these results is their remarkable simplicity, compared to the corresponding results 
in the parametrization through the Eulerian angles. 

The most familiar use of the Euler angles is of course in the description of the rotational 
motion of an orthogonal coordinate system. The instantaneous angular velocity is related to the 
time rate of change of the orientation of the rotating coordinate system by the well known equations 
[7] 

(Oi = ~~T cos <f)-\- -y- sin 6 sin </> (1 a) 

W2— "7" sin — sin 6 cos c/> (lb) 



m = —r + —rcosd, (lc) 

at at 

where oj,, co 2 , and o> 3 are the components of the angular velocity in a space-fixed coordinate sys- 
tem, and </>, 0, and i/i are the Euler angles which describe the instantaneous orientation of the 
rotating coordinate system in the spaced-fixed system. If €o(t) is regarded as known, eqs (1) are 
the differential equations governing the time development of the orientation of the rotating system. 
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With the increased use of the effective axis and angle to parametrize a rotation, it is of interest 
to examine the relation between the instantaneous angular velocity of a rotating coordinate system 
and the time development of its orientation when the orientation is specified by an effective axis 
and angle. Such an examination would give additional insight into this simple, yet rarely used 
parametrization of rotations. That examination is the purpose of this paper. 

In section 2, the method of the intrinsic vector [8, 9] is used to derive expressions for the time 
rate of change of the effective axis and angle of rotation, induced by the angular velocity. A de- 
composition is then given of the angular velocity in an orthogonal set of vectors generated by a 
unit vector along the effective axis. This decomposition is the analog of eqs (1) in the parametriza- 
tion based on the effective axis and angle. It is used to construct a new geometrical interpretation 
of the effective angle of rotation. The question of the extent to which the motion of the effective 
axis differs from the rigid rotation of an axis fixed in the rotating system, is briefly considered. It 
is shown that the motion of the effective axis fulfills the condition of rigid rotation only in the simple 
case when the rotating coordinate system rotates about a fixed axis. 




FIGURE 1. Instantaneous orientation of a coordinate system 
rotating with respect to a fixed one. 



2. Time Development of the Orientation of a Rotating Coordinate System 

We consider a pair of co-original orthogonal coordinate systems S and S' , of the same handed- 
ness, specified, respectively, by the triads of unit vectors (bi, b 2 , b 3 ) and (bi, h' 2 , bs), and take 
the point of view that S' is rotating with respect to S. The situation is shown in figure 1, in which 
co represents the instantaneous angular velocity of S'. At any instant t the triads {hi} and {b-} are 
linearly related by the orthogonal rotation dyadic A(t) according to the equation 



b;(t)=bi-A(*) = AW-b,, 



(2) 



where A is the dyadic conjugate to A. The matrix elements of A in the coordinate system S are 
given by 



Aij = bi • A • hj= hi • bj. 
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(3) 



The time derivative of A and the instantaneous angular velocity are related by 

dk ~ , 

— = coXa= (coXI) a, (4) 

in which I is the unit dyadic. Using the fact that the conjugate of co X I is -co X I, and the orthog- 
onality of A, 

A-A=A-A=I, (5) 

we may write eq (4) as 

^=-A(coXI) . 

dt (6) 

In matrix element form, this equation is 

- = v4jfc€/y m a>m, (7) 



dAtj 
dt 



where €kjm is the completely antisymmetric Levi-Civita symbol. 

The effective axis and angle of rotation for time t are extracted from A(0 by computing the 
trace of A, tr A, and constructing the intrinsic vector [9], 

Vi = €i jk Aj k . (8) 

These quantities are related to the effective angle of rotation a, and a unit vector n along the 
effective axis of rotation by the formulas 

tr A= 1 +2 cos a, (9) 

V = 2nsin a. M°) 

The handedness of the description of a with respect to n is fixed by eq (10) to agree with the common 
handedness of the coordinate systems S and S' . 

By differentiating eq (8) with respect to t, using eqs (7) and (10) to evaluate dAjkldt and dV/dt 
respectively, and making use of the identity 

€ijk €msk = Sim 8js ~ &is &jm, (11) 

we arrive at 

dn , da 

2 -j- sin a-\- 2n "77 cos a= B ■ co, (12) 

where the dyadic B is given by 

B= (1+2 cos a)\— A. (13) 

Equation (12) is the starting point for computing the equations of motion for n and a. 

We first take the scalar product of eq (12) with n. Since n is a unit vector, it is perpendicular 
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to its time derivative. Furthermore, from Euler's theorem [10] 

n • A = A • n=n. (14) 

The result of dotting eq (12) with n is therefore 

da 
2 -7— cos a= (1 + 2 cos a) co • n — co • n, 

= 2co • n cos a. 

Cancelling the common factor 2 cos a, we have 

da 

Equation (15) states the compellingly simple, and intuitively almost obvious result that the instan- 
taneous rate of change of the effective angle of rotation, is given by the component of the instan- 
taneous angular velocity along the instantaneous effective axis of rotation. 
We now insert the trigonometric representation of A [11], 

A=nn+nXlsina-nX(nXl) cos a, 

into eq (13) and use the result to evaluate B • co in eq (12). After replacing da/dt by co • n and collect- 
ing similar terms we have 

j 
2 — sinainX co sin a + n X (n X co) (1 + cos a) =0. 

This may also be written 



cfa = 1 
dt 2 



a 

nXw + nX(nXo)) cot — 



(18) 



There is no actual singularity at a = in eq (18) since for a = 0, n and co are parallel. 

Equation (15) and either of eqs (17) and (18) constitute the coupled equations of motion for 
the effective axis and angle in terms of the angular velocity. The inverse problem to that of obtaining 
these equations, is the problem of constructing the angular velocity from n, a, and their time deriva- 
tives. The construction is most readily carried out by formally solving eq (12) for co. This requires 
the dyadic B _1 which is inverse to B. To test whether the inverse exists, we compute det B, which 
is most easily done by considering the special form of A (t) when n(^) instantaneously lies along 
one of the coordinate axes, say, the b3 axis. The result of the computation is 

det B = 4 cos a (cos a+ 1). (19) 



Equation (19) shows that B _1 does not formally exist when a is 7r/2, 77, or 37r/2. In constructing 
B" 1 , we shall therefore assume that a is not equal to any of these special angles. Nevertheless, 
the final result for co will be seen to be well defined for all a. 

Using the trigonometric representation (16), and the identity 



nX(nXI) =nn- 1, 
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we may write B in the form 



^her 



B = 2 [v 2 \ - u 2 nn - uvn XI]. 



(20) 



OL 

u= sin — » 



(21a) 



a 
v= cos — • 



We then assume an expansion of B _1 of the form 

B" 1 = P\ + Qnn + Rn X I, 



(21b) 



(22) 



and compute /\ Q, and R from the requirement 

B-B" 1 = B" 1 -B= I. 
The results are 



P = 



I 



(23) 



(24a) 



<? = TT 



v L — ur 



(24b) 



R = 



2 v 



(24c) 



The expression for B l in terms of full angles 



2 



1 — cos a sin a 
I + -nnl n X I 



cos a 



1 + cos a 



(25) 



If we now multiply eq (12) from the left by B _1 , and note that n is perpendicular to its time derivative, 
and that n X I • n = 0, we arrive at 



da dn . dn , 

(ti = n — H — — sin a 4- n X — ( 1 — cos a) 
dt dt dt 



(26) 



Equation (26) is the analog of eqs (1) in the parametrization of the rotation through n and a, and gives 
the expansion of (o(t) in the instantaneous mutually orthogonal trio of vectors n, dri/dl, and 
n X dn/dt. 

An interesting geometrical interpretation of the angle a emerges from eq (26). To arrive at 
this interpretation, we first construct the unit vectors which are parallel to dn/dt and nX dn/dt 
respectively. The lengths of these vectors are equal and, from eqs (15) and (18), are given by 



dn 
dt 



n X 



dn 
dt 



( da\ 2 a 

UJ csc r 



(27) 
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The quantity [co 2 — (da/dt) 2 ] 1 ! 2 is the component of co perpendicular to n. If /3 is the angle between 
o> and n, we may write 



w2 ~(f) 2 = ft,sin/3 - 

We now define three mutually orthogonal unit vectors qi, q 2 , q3 as follows: 

dn/dt 
qi = \dnldt\' 

q 2 = n X qi, 

q 3 = n. 
In terms of these vectors and the angle ft in eq (28), we may rewrite eq (26) in the form 



(28) 



o> = co(qi sin (3 cos -~ H- q 2 sin /3 sin -~ + q3 cos ft) . 



(29a) 

(29b) 
(29c) 

(30) 



From eq (30) we see that the angle a/2 is the azimuthal angle of co with respect to the plane of qi and 
q 3 , that is, the plane of n and dn/dt. This interpretation is illustrated in figure 2, which shows the 
polar representation of o> in the coordinate system formed by (qi, q 2 , qs). 




n x dn/dt 
I dn/dt | 



FIGURE 2. Polar representation of the instantaneous angular velocity in a rotating 
coordinate system generated by a unit vector along the effective axis, and by the 
time derivative of that unit vector. 

The condition that the effective axis rotate rigidly with S' is 



dt 



co X n. 



(31) 
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From eq (18), it is clear that in general, the motion of n departs significantly from a rigid rotation. 
The condition (31) is so restrictive that the only class of motions of S' for which it is obeyed is that 
for which n is fixed in both S and S', with S' consequently rotating about the fixed direction of n. 
This can be shown directly from eq (26), but is most easily shown from Euler's theorem. Differentiat- 
ing eq (14), and making use of eq (6), we have 



dn \ dn 

— -coX n ) = — 
dt J dt 



(32) 



From eqs (31) and (32) we see that the condition that n rotate rigidly withS' leads to 



2= -0. 

dt 



(33) 



Equation (26) then reduces U 



da 
dt 



(34) 
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